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I. INTRODUCTION 



In the study of quantum gravity one encounters many technical compUcations, and it is often 
desirable to test one's ideas and tools in a simpler setting. The NonLinear Sigma Models (NLSMs) 
have striking similarities to gravity: they are nonpolynomially interacting theories, and from the 
point of view of power counting, they have exactly the same structure as gravity. On the other 
hand, they lack the complications due to gauge invariance. They are therefore a good theoretical 
laboratory where one can study various technical aspects of the renormalization of gravity without 
having to consider the complications due to gauge fixing, and with the certainty that one's results 
are not gauge artifacts. Recent work on the beta functions of gravity suggests that there might 
exist a nontrivial fixed point with finitely many UV attractive directions, making this theory 
"asymptotically safe" . This means that if one considers all the terms in the derivative expansion 
of the effective action, the corresponding (renormalized) couplings would all run towards a Fixed 
Point (FP) in the UV limit, and that only finitely many combinations of couplings would be 
relevant (attracted to the FP in the UV). Then, the requirement of tending to the FP in the UV 
would constrain the theory to lie on a finite dimensional surface, and the theory would then be 
predictive. See the original work [ll] for the definition of asymptotic safety, for reviews and [s, 4| 
for more recent results. Understanding the UV behaviour of the NLSM may shed some light on 
the analogous issue for gravity. 

Aside from this, the NLSMs also play an important role in particle physics phenomenology: they 
are used as low energy effective field theories both for strong and weak interactions. In the former 
case the scalar fields are identified with the light mesons , in the latter with the three Goldstone 
degrees of freedom of the complex Higgs doublet Q]. These effective field theories are usually 
thought to break down at some cutoff scale, of the order of the GeV in the strong case and of the 
TeV in the weak case. It is an interesting question in itself, and one that may have some relevance 
also for particle physics, whether some of these NLSM's might actually be asymptotically safe. 
Old work on the epsilon expansion and expansion suggests that a fixed point with the right 
properties may exist 0, 0, 0, More recently, the beta functions of the NLSM were recalculated 
using a two derivative truncation of an exact RG equation, and it was found in the case of the 



In the present work we begin addressing 



0{N) models that they have a nontrivial UV FP 
the issue of asymptotic safety in the NLSM taking into account also four derivative interactions. 
The beta functions of four derivative NLSM were considered before in [2l and The former 
reference uses a formalism that applies only to group-valued models; the latter uses dimensional 
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regularization and therefore cannot properly compute the running of the two derivative terms, 
which is necessary to estabhsh asymptotic safety. In this paper we extend and partly correct the 
results of these earlier works. 

This paper is organized as follows: in section 2 we discuss the models and the techniques we 
use; in section 3 we evaluate the beta functions, first for arbitrary manifolds and then for the 0{N) 
models and the chiral models; in section 4 we list their fixed points; in section 5 we close with a 
discussion. The comparison with is given in appendix A. 



II. THE THEORY 

A. Geometry and action 

In general the NLSM is a field theory whose configurations are maps from ip : X ^ Y, where 
X is a d-dimensional manifold interpreted as spacetime and Y is some n-dimensional internal 
manifold. We will always take X to be four dimensional and to have a fixed flat Euclidean metric, 
and we will call h a riemannian metric on Y. Given a map (p, one calls "vectorfield along ip" a 
rule that assigns to each point x of X a vector tangent to Y at ip{x). ^ For example, given a fixed 
vector V tangent to X at x, the image of v under the tangent map Tip is a vectorfield along ip. Its 
components are v^dfj_ip°'. Thus we can view the matrix as the components of four vectorfields 
along (p. 

The Levi-Civita connection of the metric h in TY can be used to define the covariant derivative 
of vectorfields along (p. Let Fa^y be the Christoffel symbols of h and Ra/s'^s = daFp'^'s — dpFa'^s + 
FqP eFp'^s — Fp^ ^Fa^s its Riemann tensor. The covariant derivative of a vectorfields along ip is 

v^r = 9^r+5^¥'^AV^ (1) 

A diffeomorphism / of 1" can be represented in coordinates by y' = f{y)- It maps vectorfields 
along (p to vectorfields along ip' = foip. One can check explicitly using the transformation properties 

./g _ . r' " - ^^"^ I ^^"^ ^"'^^ m 

that the covariant derivative transforms in the same way as ^ under diffeomorphisms of Y . 

We also note for future reference that the curvature of the pullback connection is the pullback 

^ The vectorfields along (p should be thought of, in geometrical terms, as sections of the pullback bundle ip*TY . 
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of the curvature of the Levi-Civita connection: 

[V;., = ^t^u^S^^ = 9^(/^"9./i2a/3^5^^ (3) 

We can now discuss the dynamics of the NLSM. Since the ordinary derivatives of (f"^ are the 
components of vectorfields along ip, the second covariant derivatives of the scalars are given by 

v^d.if'' = d^d.if'' + a^/r^%9,^^ . (4) 

Note that due to the symmetry of the Christoffel symbols V^9^(/9" = V^d^ip'^. We also define 
□(^" = D^^dfj^ip". After these preliminaries, the most general Lorentz- and parity-invariant NLSM 
with up to four derivatives has an action of the form: 

Here we defined parity to correspond to the reflection (/9"(xi, X2, X3, X4) xi, a;2, 2:3, X4). 

This is the only parity operation one can define in full generality. We will discuss below other 
"parities" that can be defined on special manifolds. At the classical level, h^'^\ h^^^ A, and T are 
fixed tensorfields on Y. They represent, in general, an infinite number of interaction terms. In the 
quantum theory these tensors will be subject to RG flow. The tensors h^'^\ /i^^^ are assumed to 
be positive definite metrics. In the present work we will always use to raise and lower indices, 
while h^"^^ is treated as any tensor. Of course nothing ultimately can depend on this convention. 
The tensor A can be assumed to be totally symmetric without loss of generality. The tensor T 
must have the following symmetry properties: 

'^I3a^5 ^a(38^ '^"/SaP • 

In ([5]) we have not considered (parity violating) terms that involve the e tensor, of the form 

where B is some four-form on Y . These could be called "Wess-Zumino-Witten terms" in a gener- 
alized sense. A proper Wess-Zumino-Witten term is one for which the four form B is not defined 
everywhere on y, but the five- form H = dB is. Then H defines a nontrivial fifth-cohomology 
class and the coefficient c has to obey a quantization condition. The original Wess-Zumino term 
corresponds to the case Y = SU{N) and H = iT{g^^dg)^. We will briefly return to these terms in 
the discussion. 
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We observe that since the field ip appears nonpolynomially in the action, it must be dimension- 
less. Then, h^'^ must have dimension of mass squared, whereas the other tensors are dimensionless. 
Later on we will find it convenient to split off a dimensionful coupling from the dimensionful tensors, 
so that all the tensors are dimensionless. 

We will be especially interested in cases in which the theory has some global symmetries. Let 
$ be a diffeomorphism of Y that leaves the tensors hP'^ ^ h^^\ A, T invariant, for example 

, , 5$^' 5$-^' 

In particular, $ is an isometry of h^^\ Then the action is invariant under the transformation 
ip ^ o ip. Such global symmetries may be discrete, or they may form a continuous group G. In 
the latter case there exist vector fields Ka on Y (with a = 1 . . . dimG) whose Lie brackets form an 
algebra isomorphic to the Lie algebra of G, and such that h^'^\ h^'^\ A, T are invariant under G: 

In particular, Ka are Killing vectors for the metric h^^^: V aKap + VjsKaa = 0. Then, the action 
([5]) is invariant under the infinitesimal transformation 5e(p" = e'^K^{ip). 

Discrete isometrics may appear in the definition of parity or time reversal. In linear 
scalar theories one can define the operation cj) i— > — For example the pions transform as 
(P7r)'*(xi, X2, X3, X4) = — 7r"(— xi, X2, ^3, a;4) under parity. In a general NLSM the transforma- 
tion ip'^ 1— > — has no intrinsic meaning. However, suppose that every point y ^ Y \s the 
:ixed point of an involutive isometry ^y. Such a manifold is said to be a symmetric space 
1^ . We can then define a new parity operation, let us call it "Parity" with capital P, by 
(P93)'^(xi, X2, X3, X4) = $0 ° ^2) 2^3) 2^4)1 where $0 is the involutive isometry of the vacuum. 

The transformation properties of the action under this new definition of parity are different than 
under the previous definition. In particular, if ^Q,^^(<I>o(y)) = ^0/37(2/)) then the j4-term will not be 
Parity-invariant. On the other hand if -60/375 ($o(?/)) = —Ba/3'ys{y), then the Wess-Zumino-Witten 
term is Parity-invariant [l^. 

B. Background field expansion 



We use the background field techniques developed in 

Q, Q, Q, Q- 



We review here some 



of the main points. Having chosen a (not necessarily constant) background ip, any other field ip 
in an open neighborhood of (p can be written y?" = ip'^ + r]°^. In principle one could work with 
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the quantum fields r/", but this is not convenient because, as differences of coordinates, they do 
not have nice transformation properties. It is therefore convenient to proceed as follows. For each 
X one can find a unique vector ^(x) tangent to (p{x) such that (p{x) is the point on the geodesic 
passing through (p{x) and tangent to ^(x), the distance between ip{x) and <^(x) being equal to 
We can thus write (p{x) = £'xp^(^).^(x), where Exp is the exponential map. The field, S,°'{x) 
is a vectorfield along (p, and its covariant derivative is defined as in ([1]). 

In principle, then, the action ip can be rewritten as S{ip) = S{(p, In practice one can compute 
the first few terms in an expansion S{(p, ^) = S^'^\(p,(,) + S^^\ip,(,) + S^'^\(p, ^) + . . . , where S^"^^ 
contains n powers of The first term is clearly S^^''{(p,^) = S{(p,0) = S((p). To compute the next 



terms we use the following formulae (whose derivation can be found in 

A bar over the derivatives and the curvatures indicates that they have to be computed with the 

background field (p. In particular for the metric g we have 

1 _ 

Inserting in ([5]), with A = 0, and keeping terms of second order in ^ we obtain 

For notational simplicity here and in the following we drop the bars over ip, V and R, but it is 
always understood that they are computed at the background field. The terms have been kept in 
the order in which they appear in ([5]), namely the first line comes from the variation of the two 
derivative term, the second and third lines come from the variation of the term containing h^^^ the 
last two lines come from the variation of the term containing T. 

C. The running effective action 

Our procedure for calculating the beta functions is a particular implementation of Wilson's 
prescription that physics at the scale k is described by an effective action F^ where all modes with 
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momenta q > k have been integrated out. We define formally an "effective average action" by 
implementing an infrared cutoff k in the functional integral over the quantum field ^. If S{ip, ^) is 
the bare action of the theory, the IR cutoff can be implemented by adding to 5 a term ASk{(p,C)j 
quadratic in which in Fourier space would have the general structure: 

ASk{'P,0 = I d\e{-q)Rkap{q')^''{q) ■ (8) 

The kernel Rkapio^)-, sometimes also called the cutoff, is chosen in such a way that the propagation 
of field modes with \q\ < k is suppressed, while field modes with \q\ > k are unaffected. There 
is a vast freedom in the choice of the cutoff ASk, and in principle physical predictions should turn 
out to be independent of this choice. One can use this freedom to simplify calculations to some 
extent. One possibility would be to write the cutoff exactly as in ([8]), with Rkapil^) = ^[f^-RfcC'?^)) 
where Rk{q^) is a scalar function of the modulus of the momentum. Note that q'^ is just the 
eigenvalue of the operator —d^ acting on the quantum field. It is more conveniente to write the 
cutoff in terms of the eigenvalues of some covariant operator, such as the Laplacian constructed 
with the background field — V^. This is the choice that was used in 



it expedient to use instead of — the full covariant fourth order operator A — • 



111 ]. In this paper we will find 
A5fc(^,0 = \j d'xeh^^l{^)Rk{A)e . (9) 



Because A depends only on the background field, and not on the quantum fields, this cutoff is still 
quadratic in the quantum fields, as required. 

Having modified the propagator of the theory, we define a generating functional 
depending on the background (p and on a source field j coupled linearly to the quantum field ^, by 

Wk{v,j) = - logy (dr)exp {-S{^,i) - ASki^X) - j J^r) (10) 
Then we define a modified /e-dependent Legendre transform 

rfc(^,0 = Wk{v,j) - j ja^ - ASk{v,0 , 

where ASk has been subtracted. The "classical fields" are denoted again for notational 
simplicity. The functional F^ reduces for /c — > to the usual background field effective action 
T{ip,^), the generating functional of one-particle irreducible Green functions of S,. 
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D. The one loop beta functional 



At one loop one can evaluate the functional F^: 



Note that AS^ has cancelled out. The only remaining dependence on k is in R^, so 



k- 



dT 



k 



1^ / 6^S „ \ , dRk 



dk 2 

The r.h.s. can be regarded as the one loop beta functional of the theory. The individual beta 
functions can be read off by isolating the coefficients of various operators. Of course one could 
derive the one loop beta functions in other, more traditional ways. We prefer this route because 
it has a few advantages. First, due to the rapid fall off of the function kdkRk, the beta functional 
is itself finite and one does not actually need to introduce any ultraviolet regularization. So, even 
though the derivation of the equation from a functional integral was formal, because the functional 
integral is itself ill defined, the functional RG equation is itself perfectly well defined. A second 
important point is that the RG improvement of this equation, where one replaces the bare action 



20|. So although in the present work we 



S by Ffc in the r.h.s., is actually an exact equation 
shall restrict ourselves to the one loop approximation, the formalism is ready for the calculation of 
the beta functions based on a truncation of the exact RG equation, which amount to resumming 
infinitely many orders of perturbation theory. A final, important point is that experience with 
other systems shows that this procedure gives exactly the same results as any other procedure for 
the universal (scheme-independent) one loop beta functions. We will see in section HID that, to 
the extent that a comparison is possible, this expectation will be confirmed also in this case. 



E. Global symmetries 

If there are any symmetries, one can define the RG flow so as to preserve them. To see this, 
let ^ be an internal symmetry, as in section IIA. Since it is an isometry of h^^\ it also leaves the 
connection invariant, so it maps the geodesic through y tangent to ^ to the geodesic through $(y) 
tangent to r$(0 UM' 
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We call ip' = ^ o (f and = r<I>(^) the transform of (p and under Then if' = ^{Exp^^) = 
Exp^{ip){T^{C)) = Exp^i^'. There follows that 



i.e. the background field action S is G-invariant provided both background and quantum field are 
transformed. The operator A is covariant, so A'(^') = T<I>(A(^)) or abstractly A' = T<I> o (A) o 
T$~^, so also the cutoff term @ is invariant: 



One can formally choose the measure in the functional integral (jlOp to be invariant under Since 
both measure and integrand are invariant, the effective action F^. will also be invariant, for all k. 

Somewhat less formally, one can arrive at the same conclusion as follows: observe that the 
cutoff as defined in ([9]) is a suppression term that depends on the eigenvalue of the operator A 
on the normal modes of the field. From the transformation properties of A one sees that if ^ is 
an eigenvector of A with eigenvalue A, then ^' is an eigenvector of A' with the same eigenvalue. 
Therefore the spectrum of A is invariant. Equation (jl2p gives the (one loop) scale variation of 
Ffc((/?) as a sum of terms, each term being a fixed function evaluated on an eigenvalues of A. Since 
the eigenvalues are invariant, the sum is also invariant, so it follows that 9fFfc((/7) is invariant. 
This implies that if the starting action TkQ^ip) is invariant, also the action at any other k is. This 
argument is mathematically more meaningful, because unlike the one based on the path integral, 
it involves only statements about finite expressions. 

The previous argument can be applied both to discrete and continuous symmetries. For example 
in the case of discrete symmetries, it implies that the fiow preserves Parity. If the A term violates 
Parity, it must be set to zero in order to have a Parity invariant theory. The flow will preserve 
this property, so the beta function of A will be zero. In other words the condition A = will be 
"protected by Parity". We will see this in an explicit calculation in section Illb. 

III. EVALUATION OF BETA FUNCTIONS 

The one loop RG flow equation (jl2p can be approximated by resorting to a truncation, which 
means keeping only a flnite number of terms in F^, inserting this ansatz in the flow equation and 
deriving from it the beta functions of the couplings that enter in the ansatz. The best way of 
truncating F^ is to do so consistently with a derivative expansion, i.e. to keep all the terms with 



(13) 



A5fc((/p',0 = A5,.(v9,e) 



(14) 
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a given number of derivatives. In this paper we will approximate Tk by a functional of the form 
(0), where the tensors h^'^\ h^^^ and T are fc-dependent, and ^4 = 0. In general this is still a 
functional flow, because the tensors actually contain infinitely many couplings. We will be able to 
say more in the case when a global symmetry restricts the possible form of these tensors, so that 
only finitely many couplings remain. In this paper we will explicitly compute the beta functions 
in the case when y is a sphere or a special unitary group. Since these are symmetric spaces, it will 
be consistent to neglect the A terms altogether. 

A. The inverse propagator 

Integrating by parts one can rewrite ([7]) in the form S^'^\(p,(,) = ^(^,A^), where the inner 
product of vectorfields along if is (^, C) = / d"^^ ^c^^"C^ ^-iid A is a self-adjoint operator of the 
form: 

A^p = h^^^a' + e^^v^y, + c^pV^ + . (15) 

Self-adjointness means that {^,A.C,) = (A^,(^) and implies the properties: 

K'p = , (16) 
C^p = -C^^ + V.e;,: + V^BJ^ , (17) 
V^fi = Vp^ + V,C^„ + VuV.Bj^ . (18) 

In addition by commuting derivatives we can arrange the operator so that = B'^^^. In order to 
arrive at the operator A we proceed in two steps. First we put all the derivatives of ([7]) on one of 
the ^'s, so that = ^(C, A^), where A is of the form ((T5|), with 

= ^''"(-^S + 29p^^aV'(i?a7/35 - TaMs)) - ^d^V^d,ip'T^-,ps , (19) 

We do not display the form of C^^, since it does not contribute to the expressions we want to 
calculate, as will become clear in due course. This operator A is not self-adjoint, and we define 
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A = ^(A + A'''). Its coefficients are 



2 V /3" 



a/3 



13a 



(21) 



Note tliat the last two terms in and are total derivatives, and will not contribute to our 
final formulae. Finally we can symmetrize i3^^ in ^, v at the cost of generating a commutator term 
that contributes to T)ap. The final form of the operator A is (jlSp . with 



V. 



a/3 = 2 + ^/3"j - dpip^dP^'d^ip'd^ip''iT^^,^Rsr,(3^ + Tp^e^Rsr,J) + TD , (22) 



where TD stands for "total derivatives" . Again we omit to give C ^g, because it does not contribute 



to the beta functions. These formulae agree with (3.17-21) in [ij], except for a factor 2 in the 
coefficient of the first term containing T^p^^ in equation (|19p. 



B. Beta functionals 



We begin by discussing the general case of the action ^ with arbitrary 

/i(4) and T, and 

^ = 0. We evaluate the trace in (jl2p by heat kernel methods. The advantage of this procedure is 
that pieces of the calculation are readily available in the literature. Given a differential operator 
A of order p, and some function W , we have 

1 



TrT^(A) 



(4^) 



Q±{W)Bo{^) + Q2{W)B2{^) + Qo(VF)B4(A) 



The heat kernel coefficients are defined by the asymptotic expansion 



Tr(e 



(47r)2 



(23) 



(24) 



with Bn = f d'^xtvbn', bn are matrices with indices a, f3 and tr denotes the trace over such indices. 



For a fourth order operator of the form ()15p . they can be found in [24]. The quantities QniW) 
in (USD are given by Qn{W) = ^ dzz^'-^Wiz) for n > and Qo{W) = W{0). We do 
not need any higher coefficients. In order to be able to evaluate the integrals in closed form we 
choose the "optimized" cutoff function Rk{z) = (A;^ — z)6{k'^ — z) 23|]. The scale derivative of 



the cutoff is = Ak'^6{k'^ — z), and the modified inverse propagator Pk{z) = z + Rk{z) is 
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equal to k'^ for z < k^. Then the function to be traced in the ERGE is just a step function: 
W{z) = |-^fc^r = 26(1 — z/k'^), and the integrals are very simple: 

4 



Qi = 2k\ Qi=^k' 

2 J-K 



Qo 



(25) 



The first term in (j23p is field independent and will be omitted. Putting together the remaining 
pieces: 



k- 



^ I d'x (^-k^Bi + lnf^n$: + ^^Bf^BZ + ^SapS"' - (26) 



dk (4t 

where il. is defined as in ([3]) and B = B^. The first term comes from B2, the others from B4. One 



finds 



1^" 

4 ° 



1 



24 A''^ po 4g "/^ 



(27) 
(28) 



S<5 



(29) 



/ X 



(30) 



From here one can read off the beta functionals of /i^^^ , yl, T as the coefficients of terms containing 
two, three and four powers of d^ip", respectively. We do not give these general formulae, but just 
make some observations. The only term proportional to □(^"□(^^ is contained in —"D^, so the beta 
functional of h^^^ is easily obtained: 



k4rh 



(4) 



dk 87r2 

This is very similar to the result for the two derivative truncation. In order to 



(31) 



compare results 



obtained with the same type of cutoff, we should repeat the calculation of [Ul] using a cutoff 
constructed with the full inverse propagator Aap = —h^^^V"^ — d^ip'' df^ip^ Ra-fi35 ■ This is a cutoff of 
type III in the terminology used in [3] . In this case the general beta function of the metric is 



k4rh^'^ 



1 Rk] 1 2 

dk"»^ - Ji^^' i 7^ ) - 8^^ ^'^^ 



(32) 
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(2) 

where R denotes now the curvature of h^p- As a side remark, this little calculation is also useful 
to test the scheme dependence of the results: with the type I cutoff used in the result was 



^^^S - (4^^2 Rap - Y^k'Rap , (33) 

which differs by a factor 2. 

Another fact that follows from ()30p is that the beta function of A (coming from the coefficient of 
\Jnp'^d^(p^ d^ip'^) is proportional to covariant derivatives of the Ricci tensor. For symmetric spaces 
the covariant derivative of the curvature vanishes and therefore on such spaces it is consistent to 
set A = Q. This confirms the general statement made in section HE. The particular models that 
we shall consider in the following are symmetric spaces. 



C. The spherical models 

We now consider the class of models for which the target space Y is the sphere S^. Such 
models are often called the 0{N) models, with N = n + 1, because they have global symmetry 
0{N). There is only one 0(n + l)-invariant nonvanishing rank two tensor on the sphere, there is 
no invariant rank three tensor and there are only two invariant rank four tensors with the desired 
index symmetries, up to overall constant factors. If we regard 5" as embedded in R"+^, we call 
haf3 the metric of the sphere of unit radius. Its Riemann and Ricci tensors are given by 

Raf3^5 = ha^hjss — h^ihp^ ; Rap = {n — l)hai3 ; R = n{n — 1) . 

Therefore both /i^^) ^nd h^^^ must be proportional to h, and T is a combination of /I's: 

Here has mass dimension 2, while A, ii, £2 are dimensionless 2. It is convenient to regard j as 
the overall factor of the fourth order terms; then we define the ratios between the three coefficients 
of the four-derivative terms as /i = A^i and /2 = Xi2- For the reader's convenience we rewrite the 
action of the 5" models: 



(34) 



the names ii and I2 are used commonly in chiral perturbation theory 0|. 
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One then finds the following beta functions: 

= -^A^ (35) 



87r2 

Ph = 4^ ((^ + + 20/2/1 + 4/1 + 6(n + 3)/i + 24/2 + 8) (36) 



Pf. = ^ (^/f + ^/i/2 + ^/l - (n + 3)/, - (3n + 1)/. + n - I) (37) 

~4 \ ~2 

Pg2 = 2f + r^ ((5 + n)/i + (2 + 4n)/2 + 4(1 - n)) - ^ ((5 + n)/i + (2 + 4n)/2 + 2(1 - n))(38) 

Equations (j36p and (j37p differ in a significant way from equation (5.11) in [l^. This is due to the 
already mentioned factor 2 in a term in pop . Unfortunately, this changes completely the picture 
of the fixed points. 

It will be instructive to compare the results of this four derivative truncation with those of the 
simpler two derivative truncation discussed in llj. If we specialize ()32p to y = /S", it gives 



whereas from (138p . setting for simplicity ii = £2 = and in the limit A — > one gets 

The difference is just a factor 2, which is within the range of variation due to the scheme dependence. 
It is quite remarkable that the beta function is so similar in spite of the very different dynamics. 
We shall see in section IVA that this fact is quite general. 



D. The chiral models 



Next we consider the case where Y is the group SU{N). In this case it is customary to denote 
U{x) the matrix (in the fundamental representation) that corresponds to the coordinates We 
demand that the theory be invariant under left and right multiplications U{x) ^ gJ^^U{x)gFi, 
forming the group SU{N)l x SU{N)ji ("chiral symmetry"). Further we demand that the theory 
be invariant under the discrete symmetries U{x) 1— > [/"^(x), which corresponds physically to charge 
conjugation, to the simple parity xi ^ —xi, to the involutive isometry : U ^ and hence to 
Parity [/(xi, X2, X3, 0:4) 1— > U~^{—xi,X2,X3,Xi). More details on the translation between the tensor 
and the matrix formalism are given in Appendix A. 

Let Ca be a basis of the Lie algebra, with a = 1 ... — 1. We denote Tq the corresponding 
matrices in the fundamental representation; they are a set of hermitian, traceless N x N matrices. 
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We fix the normahzation of the basis by the equation 

TaTb = -^^ah + ^ {dabc + ifabc) Tc ■ (41) 

(In the case of SU{3) these matrices are one half the Gell-Mann A matrices.) 

A tensor on SU{N) which is invariant under SU{N)l x SU{N)ii is said to be "biinvariant" . 
There is a one to one correspondence between biinvariant tensors on SU (N) and ^dd-invariant 
tensors in the Lie algebra of SU (N), where Ad is the adjoint representation. Given an j4(i-invariant 
tensor tah...'^'^'" on the algebra, the corresponding biinvariant tensorfield on the group is 

J. 7(5... J. cd...rarb TlT^ 

^ap... — ^ab... ^a^(3 ■ ■ ■ ^c^d- • • 

where are the components of the left-invariant Maurer Cartan form L = U~^dU = L'^dy"^ {—iTa) 
and are the components of the left-invariant vectorfields on SU {N). The matrix is the inverse 
of L^. (In this construction we could use equivalently right-invariant objects.) 

Up to rescalings, there is a unique j4(i-invariant inner product in the Lie algebra, which we 
choose as hab = 2Tr TaT^ = 5ab ^- Then the corresponding biinvariant metric is 

V = KL^Kb , (42) 

so that the left-invariant vectorfields La can also be regarded as a vierbein. The Riemann and 
Ricci tensors and the Ricci scalar of h are given by 

RaP^S = -^L'^L^fjL'^Lj fab'^ feed Rafi = -^NhajS ] R = -N{N'^ - 1) . (43) 

As with the sphere, we define h^^^ = -^hap, h^^^ = jhap- The tensors dabc and fabc are a totally 
symmetric and a totally antisymmetric ^d-invariant three tensor in the algebra. In principle chiral 
invariance would permit a term in the action with A^p-y = L'^L^^L^dabc', however using L^($o(y)) = 
Raiy), La{y)Ra{y) = Ad{g{y)fa and the ^d-invariance of dabc, one sees that Aapyi^oiv)) = 
Aaf3y{y), so this term violates Parity. 

For T we have the following Ad-invariant four-tensors in the algebra with the correct symmetries: 

^abcd 2 ^^"•^^^'^ ~^ ^ad^bc) j ^abcd ~ ^ab^cd j '^abcd ~ 2 ^fo,cefbd ~l~ fadefbc ) j 

T^ld = \ idaccdbd' + dadcdbcl T^tld = dabedcd" • (44) 

^ Here the matrices are in the fundamental representation. The Cartan-Killing form just differs by a constant: 

Bab = Tl{Ad{Ta)Ad(Tt)) = NSab- 
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They are not all independent, however. The identity (2.10) of 2j] implies that 



j^T^^^ - ^T^^^ + T'-^^ + ^^^^ - T^^^ = , 

so that r(5) can be eliminated. In the case = 3 the identity (2.23) of 
preceding relation, further implies 



(45) 



24| . together with the 



(46) 



so that we can also eliminate T^^^ . Finally in the case N = 2 the tensor date is identically zero, so 
we can keep only T^^^ and T^^^ as independent combinations, and use T^'^^ = T*-^-* — T^^\ 
The action of the generic SU (N) models can then be written in the form: 

4 

" ' ^ 0/375 



1 



i=l 



(47) 



and the sum stops at i = 3 and i = 2 for = 3 and N = 2 respectively. As in (|34p . it will be 
convenient to use instead of the couplings ij the combinations fi = Xii . 



Making repeated use of traces given in 



25I ] one finds the following beta functions: 



327r2 
A 

7687r2Ar2 



(48) 



16A^2(Af2 + 20) ff + 64N^f^ + ISOA^VI + 4(149 A^^ - 1280)/| 
-320Ar2/i/2 - 32Ar3/i/3 + 32Ar(A^2 ^ ^^j^f^ ^ l28Nf2h - l^ON^hh 
-24Af3/i - lOSA^Vs + 36Af2/4 + 9A^^ 



(49) 



A_ 

7687r2Af2 L' 



8Af^(A^^ + 14)/f + 32Ar^(6Af^ + l)f^ + GON^f^ + 4(7A^^ + 656)/| 



+32Ar2(3Ar2 + 14)/i/2 + 80Ar3/i/3 + 16Ar(7Ar2 - 44)/i/4 + 288Ar3/2/3 



(50) 



+32Ar(15Ar2 - 64)/2/4 + 120Ar2/3/4 - 2AN\fi + 3/2) - 36Ar2(/3 + f^) + 3N^ 



Ph 



Ph 



Po 



A 



15367r2A^ 



52ArV| + 12(23Af^ - 320)/| + JdSNfifs + 2hm hh + ^^Nfih + l28Nf2U 



+24(llAr2 - 64)/3/4 - 192Af(/i + h) - 60N\h + h) + 384/4 + A^^ 



A 



15367r2A^ 



GON^fl + 4(87Ar^ - 1728)/| + 1536Ar/i/4 + 768Ar/2/4 + 21QN'' hU 



-36N\h + h)+SN'' 



2f + {N{N^ + 4)/i + 2N{2N^ - l)/2 + 3Ar2/3 + - 4)/4 - A^') 

~Tm^ + + 2Af(2Ar2 - l)/2 + 3Ar2/3 + 5{N^ - A)U - /2) 



(51) 



(52) 



(53) 



In appendix A we establish the dictionary between our notation and that used in [Ij]. When the 
beta functions are compared, we find perfect agreement, except for one small difference: the very 
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Q, i.e. 



and would have the 



last term in the first line of f5g2 would be N /2 according to 
same coefficients. This is the same difference that we observed between (j32p (type III cutoff) and 
(j33|) (type I cutoff Vso, effectively the calculation in [3l is equivalent to a type I cutoff. Given that 



the calculation in [1^ was done using completely different techniques, this agreement confirms that 
the one loop beta functions of the dimensionless couplings (which in a calculation of the effective 
action would correspond to logarithmic divergences) is scheme independent. 

The cases = 3 and N = 2 have to be treated separately, because in these cases only three, 
respectively two, of the couplings fi are independent. In the case = 3 one can eliminate in 
favor of the other three couplings. Then using (j46p one can obtain the beta functions of /i, /2 and 
/a from the ones given above by 

A 



N=3 
N=3 

A 



Ar=3,/4=0 



7687r2 



[464/2 + 64/1 + 180/1 + 320/1/2 - 96/1/3 + 72/1 - IO8/3 + 9] 



15367r2 

■'^ N=3 



N=3,f4=0 
2 I ion^2 



368/i^ + 3520/2" + 180/1 + 2624/1/2 + 48O/1/3 + 1728/2/3 - 144/1 - 432/2 - IO8/3 + 9 



1 



A 



[2/1 + I6/1/3 + 8/2/3 - 4/1 - 4/2 - 3/3 



Ar=3,/4=0 327r2 

In the case N = 2 we can set = 0, because T^^^ = identically, and we can eliminate f^. 
One can obtain the beta functions of /i, /2 from the ones given above by 

A 



N=2 



N=2 



Pfi - Ph 

Ph + Pfs 



Ar=2,/3=0,/4=0 967r2 

A 



1927r^ 



[48/f + 8/1 + 40/1/2 + I8/1 + 12/2 + 1] 
[36/2 ^ 200/1 + 2O8/1/2 - 36/i - 6O/2 + 1] 



N=2,h=0j4=0 

The latter result can be used to check also our beta functions for the spherical sigma model. In 
fact there is exactly one manifold which is simultaneously a sphere and a special unitary group: 
it is SU{2) = S^. Thus the beta functions should agree in this case. Before comparing, a little 
point needs to be addressed. In section IIIC we chose the metric ha/s to be that of a sphere of 
unit radius. In this section we have fixed the metric by the conditions (j4ip . (j42p . It turns out that 
in the case N = 2 this normalization corresponds to a sphere of radius two. This can be seen for 
example from equation (j43p . specialized to = 2, with fabc = ^abc- In order to compare the beta 
functions of 5'^ with those for SU{2) we therefore have to redefine A — > A/4, /i — > 4/i, /2 4/2, 
— > 5^/4. With these redefinitions, the beta functions do indeed agree. 
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IV. FIXED POINTS 



A. The spherical models 

We now discuss solutions of the RG flow equations. The beta function of A depends only on A 
and the solution is 

^(') - TTW^ • 

where Aq = A(to)- We assume Ao > 0, thus A is asymptotically free. The beta functions of fi 
and /2 do not depend on g, so their flow can be studied independently. Here we do not discuss 
general solutions but merely look for fixed points. The overall factor A in these beta functions can 
be eliminated by a simple redefinition t = t{i) of the parameter along the RG trajectories: 

-- = T#- (55) 

Since i is a monotonic function of t, the FPs for fi and /2 are the zeroes of the modified beta 
functions 

They are just polynomials in /i and f2- The model has no real FP for n = 2, but there are FPs 
for all n > 2. For n = 3, . . . 8 they are given in the fifth and sixth column in Table I. One can then 
insert the FP values of /i and /2 in (3g2 and look for FP of g"^. In each case there are two solutions, 
one at g^ = 0, the other at some nonzero value. These solutions are reported in the fourth column 
in Table I, for n = 3, . . . 8. The first solution describes the Gaussian FP (GFP), where all the 
couplings g'^, A, 1/^2 are zero, the others non Gaussian FP's (NFP) where g^ has finite limits 
instead. Each FP can be approached only from specific directions in the space parametrized by A, 
ii, £2, i-e. the ratios /i and /2 take specific values. For each NFP these values are unique, while 
for the GFP there may be several possible values: two if n = 3, 4, 5 and four if n = 6, 7, 8. 

When one considers the linearized flow around any of the GFPs, one finds as expected that the 
critical exponents, defined as minus the eigenvalues of the matrix are -2,0,0,0, corresponding to 
the canonical dimensions of the couplings. The critical exponents at the NGP are instead 2,0,0,0. 
Thus the dimensionless couplings are marginal, and of the two FPs, the trivial one is IR attractive 
and the nontrivial one UV attractive for g. For A it is clear that the FP is UV attractive (if we had 
chosen A < it would be IR attractive). In order to establish the attractive or repulsive character 
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of /i and /2, one can look at the linearized flow in the variable t, which is described by the 2x2 
matrix 




We define the "critical exponents" 61^2 to be minus the eigenvalues of this matrix. They are 
reported in the last two columns of table I, for n = 3, . . . 8. It is important to realize that even for 
the GFP the eigenvectors of the stability matrix are not the operators that appear in the action 
but mixings thereof. We do not report the eigenvectors here. 

Beyond the values given in Table I, we have checked numerically the existence of the FP up 
to n = 200. For large n one can study the theory analytically, to some extent. There are four 
FPs for the system of the /j's, which are: /i = 0, /2 = 1 with critical exponents 9i = 6, O2 = 12, 
/i = 0, /2 = 1/2 with critical exponents 0i = 6, O2 = —12, /i = —6, /2 = 5/2 with critical 
exponents 9i = —6, 62 = 12, fi = —6, /2 = 2 with critical exponents Oi = —6, 62 = —12. The 
numerical values at finite n do indeed tend towards these limits for growing n. 



B. The chiral models 



The chiral model with A'^ = 2 is equivalent to the spherical model with n = 3 (up to the 
redefinition of the couplings mentioned in the end of section HID) so we need not discuss this 
case further. For ease of comparison we just report the properties of its nontrivial FPs in the 
parametrization we used for the chiral models: 

NFPl : fu = -0.173 ; /2* = 0.113 ; g = 13.25 
NFP2 : fu = -0.261 ; /2* = 0.154 ; g = 12.78 

The critical exponents don't depend on the definition of the couplings and therefore are the same 
as in Table I; they do however depend on the choice of RG parameter and they differ from those 
given in IJ] by a factor 47r^, which is due to the definition of the parameter x there. 



In the case = 3 the system of the /j's has two FPs at 

FPl : fu = -0.154 ; /2* = 0.050 ; /g* = 0.085 ; 
FP2 : fu = -0.108 ; /2, = 0.043 ; /g, = 0.061 . 

The attractivity properties in the space spanned by the /j's is given, as in the spherical case, by 
studying the modified flow with parameter t. The critical exponents at FPl are: 0.0303 with eigen- 
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n 


9* 


FP 


9* 


/i* 


/2* 


Oi 


02 


3 


8.886 


NFPl 


6.626 


-0.693 


0.453 


0.094 


-0.0121 


3 




NFP2 


6.390 


-1.042 


0.615 


0.103 


0.0119 


3 




GFPl 





-0.693 


0.453 


0.094 


-0.0121 


3 




GFP2 





-1.042 


0.615 


0.103 


0.0119 


4 


7.255 


NFPl 


5.877 


-0.479 


0.398 


0.105 


-0.0412 


4 




NFP2 


5.442 


-1.555 


0.852 


0.132 


0.0392 


4 




GFPl 





-0.479 


0.398 


0.105 


-0.0412 


4 




GFP2 





-1.555 


0.852 


0.132 


0.0392 


5 


6.283 


NFPl 


5.310 


-0.400 


0.400 


0.118 


-0.0608 


5 




NFP2 


4.924 


-1.875 


0.988 


0.154 


0.0567 


5 




GFPl 





-0.400 


0.400 


0.118 


-0.0608 


5 




GFP2 





-1.875 


0.988 


0.154 


0.0567 


6 


5.620 


NFPl 


4.883 


-0.350 


0.408 


0.131 


-0.0780 


6 




NFP2 


4.577 


-2.131 


1.091 


0.171 


-0.0717 


6 




GFPl 





-0.350 


0.408 


0.131 


-0.0780 


6 




GFP2 





-2.131 


1.091 


0.171 


0.0717 


6 




GFP3 





-0.814 


1.369 


-0.161 


-0.0539 


6 




GFP4 





-2.363 


2.091 


-0.164 


-0.0617 


7 


5.130 


NFPl 


4.548 


-0.314 


0.417 


0.143 


-0.0939 


7 




NFP2 


4.322 


-2.347 


1.175 


0.185 


0.0851 


7 




GFPl 





-0.314 


0.417 


0.143 


-0.0939 


7 




GFP2 





-2.347 


1.175 


0.185 


0.0851 


7 




GFP3 





-2.790 


2.130 


-0.181 


-0.0647 


7 




GFP4 





-0.598 


1.241 


-0.174 


-0.0716 


8 


4.750 


NFPl 


4.274 


-0.286 


0.424 


0.156 


-0.1092 


8 




NFP2 


4.125 


-2.535 


1.247 


0.197 


0.0976 


8 




GFPl 





-0.286 


0.424 


0.156 


-0.1092 


8 




GFP2 





-2.535 


1.247 


0.197 


0.0976 


8 




GFP3 





-2.790 


2.131 


-0.180 


0.1023 


8 




GFP4 





-0.598 


1.247 


-0.187 


-0.0872 



TABLE I: Gaussian and non-Gaussian fixed points of the S"" model at one loop. The first column gives the 
dimension ?i. The second column gives the position of the NGFP in the two-derivative truncation, using a 
type III cutoff. The rest of the table refers to the four-derivative truncation, also using a type III cutoff. 
The third column gives the name of the FP. Columns 4,5,6 give the position of the NGFP, columns 7,8 the 
critical exponents, as defined in the text. The coupling A, not listed, goes to zero and is marginal in this 
approximation. 
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vector (0.411, 0.630, 0.658); 0.0123 with eigenvector (0.515, -0.570, 0.640); 0.00289 with eigenvec- 
tor (0.869, -0.148, -0.473), whereas at FP2 they are: 0.0280 with eigenvector (0.366, 0.618, 0.695); 
0.0108 with eigenvector (0.513, -0.575, 0.638) and -0.00293 with eigenvector (0.887, -0.125, -0.445). 
Therefore FPl is attractive in ah three directions, while FP2 is attractive in two directions. For 
each of these two FP's, the beta function of g has two FP's: the trivial FP, which has always 
critical exponents -2, and a nontrivial FP, which is located aX g = 11.17 for NFPl or 11.50 for 
NFP2, and having critical exponent 2 in both cases. 

We have found no FP's for iV > 3: the system of equations = for i = 1,2,3,4 only 
has complex solutions. To cover all of theory space we have checked this statement also in the 
parametrization of the (.i and in the parametrization of Ui = l/^j. This is true also in the large N 
limit. If we keep only the leading terms (of order A^^ for /i and /2 and of order N for /s and 74), 
again the resulting polynomials do not have any real zero. 



V. DISCUSSION 



We have calculated the one loop beta functionals of the NLSM with values in any manifold, in 
the presence of a very general class of four derivative terms. We have then specialized our results 
to two infinite families of models: the 0{N) models, with values in spheres, and the chiral models 
with values in the groups SU{N). Such calculations had been done before, but since the results are 
rather complicated, it is useful to have independent verifications. Our approach is calculationally 
very similar to Q, but after correcting some small errors at the general level, we find that the 
FP structure of the 0{N) models is completely different from their findings. On the other hand 



121 for what concerns the dimensionless 



our results for the chiral models agree completely with 
couplings, even though the calculation was done using very different techniques. Since SU (2) = S^, 
this provides a check also for our results for the spheres. 

Since our aim is to establish asymptotic safety, or lack thereof, it is important for us to have 
also the beta functions of the dimensionful coupling g, which in the chiral models is the inverse of 
the pion decay constant. This had not been considered at all in [3], but it had been calculated in 
12 1 for the chiral models. Again we have agreement with the result of [3], up to a single factor 2 



in one term; as discussed before, since this beta function is scheme dependent, we believe that this 
is not an error on either side, but the result of the different way in which the calculation was done. 
This difference results in a shift of the FP value of g; for example in the case of SU{2) one would 
find g = 19.88 instead of 13.25 for NFPl and 18.39 instead of 12.78 for NFP2. Such variations by 
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a factor of order 2 are to be expected. 

One of the motivations of this work was to use the NLSM as a toy model for gravity. From this 
point of view we have a perfect correspondence of results. If we use the 1 j-p^ propagator that comes 
from the two derivative term, both theories are perturbatively unitary but nonrenormalizable; if on 
the other hand we use the propagator that comes from the four derivative terms both theories 
are renormalizable (see [26[ for gravity and 



case it had also been established (see 



28 



29|, 



2711 for the NLSM) but contain ghosts. In the latter 



3d, 



3l| for gravity and 



Q,Q] 



for the NLSM) that 



the four derivative terms, whose couplings are dimensionless, are asymptotically free. Actually 
the analogy works even in greater detail. The coefficient of the square of the Weyl tensor (for 
gravity) and the square of Dc^" (for the NLSM) have at one loop a beta function that is constant. 
These coefficients diverge logarithmically in the UV, so their inverses, which are the perturbative 
couplings, are asymptotically free. The coefficients of the other four derivative terms have more 
complicated beta functions, but overall there is asymptotic freedom, provided the Gaussian FP is 
approached from some special direction. There have been many attempts to avoid the effects of 
the ghosts, see 



28 



321 ] for gravity and 



12l | for the NLSM. In any case, the existence of the ghosts is 
only established at tree level. Whether they exist in the full quantum theory is a deep dynamical 
question whose answer is not known. 

All these "old" works on higher derivatives theories concentrated on the behavior of the couplings 
that multiply the four derivative terms; much less attention, if any, was payed to the coefficient of 
the two derivative term, which has dimension of square of a mass: the inverse of Newton's constant 
in gravity and the square of the pion decay constant in the chiral NLSM. In several papers this 
issue was ignored, or incorrect results were given, because of the use of dimensional regularization. 
The correct RG flow of these couplings is quadratic in A;, and is best seen when a momentum cutoff 
is used. In [33] this point was made for gravity and it was shown that when this quadratic running 



is taken into account the beta function for Newton's constant (and for the cosmological constant) 
has, in addition to the Gaussian one, also a nontrivial FP. In this paper we have found that the 
same is true in a large class of NLSM. This is crucial for asymptotic safety. 

It is somewhat gratifying to see that the FP does not always exist for all NLSM: in particular 
we have seen that within the one loop approximation, adding the higher derivative terms destroys 
the FP that is present in the two derivative truncation for the sphere 5^ and for the chiral models 
with > 3. If there was any doubt, this shows that the existence of the FP is not "built into the 
formalism" but is a genuine property of the theory. This is somewhat analogous to the situation 
when one adds minimally coupled matter fields to gravity 34l ]. 
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The next step will be to replace the one loop functional RG equation (jl2p by its exact coun- 
terpart, which only differs in the replacement of the bare action S by in the r.h.s. [2^. There 



that needed further 



are at least two good reasons to do this calculation. One of the points of 
clarification was the value of the lowest critical exponent. In the two derivative truncation at one 
loop it was always 2 at the nontrivial FP. Thus the critical exponent that governs the rate at 
which the correlation length diverges was given by 

1 _ 1 _ 1 

at * 

which is the value of mean field theory. Using the "exact" RG truncated at two derivatives gave 
u = 3/8 for the 0{N) models, independent of N . One would like to understand what effect the 
higher derivative terms have on this exponent. Since here we restricted ourselves to one loop, 
we found again v = 1/2, so the calculations of this paper are of no use in this respect. Another 
motivation comes from recent calculations in higher derivative gravity [J] that go beyond one loop 
and find that the theory is not asymptotically free, but rather all couplings reach nonzero values 
at the UV FP. It would be interesting to see similar behaviour in (some) NLSM. 

Concerning possible direct phenomenological applications of the NLSM, regarded as an effective 
field theory, it is interesting to ask what relation, if any, the UV properties of the NLSM may have 
to the properties of the underlying fundamental theory. Regarding the chiral NLSM as the low 
energy approximations of a QCD-like theory, one may note that there is rough agreement between 
the range of existence of the NLSM FP and the "conformal window" for the existence of an IR 
FP in the case when the quarks are in the adjoint or in the symmetric tensor representation 35l |. 



One could get a better understanding of this issue if the beta functions of the NLSM depended 
on the number of "colors" of the underlying theory, which in the effective theory are reflected in 



the coefficient of the Wess-Zumino-Witten term 
Zumino-Witten term is zero 



13 



15( |. The one loop beta function of the Wess- 



3g]; this is consistent with the quantization of the coefficient c. 
Unfortunately the beta functions of the remaining couplings are completely independent of this 
coefficient, so the low energy theory seems to be insensitive to this parameter. 



Another possible application is to electroweak chiral perturbation theory jg]. If the NLSM 
turned out to be asymptotically safe in the presence of gauge fields and fermions, then one may 
envisage a higgsless standard model up to very high energies. This will also require a separate 
investigation. A related application of asymptotic safety to the standard model has been discussed 
recently in [s^]- 

To summarize we believe that the NLSM are interesting theoretical laboratories in which one 
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may test various theoretical ideas, and they have also important phenomenological applications. 
The question whether some NLSM could be asymptotically safe seems to us to be a particularly 
important one, and to deserve more attention. 
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VI. APPENDIX A 



In 12] the action for the chiral SU{N) model is written in the form: 



■^csTriL^L^'UL'' + L^L^L'^L") - C4T¥(L^L'^)Tr(L^L^) - C5Tr(L^,L'^)TV(L^L'') . (56) 



where = {/ ^df^U. We want to translate this action into the form (I47p . Deriving the equation 
= d^ip°'L%{-iTa) we obtain 

The antisymmetric part of this equation is 

d^Ly — dyL^ = — [L^jLu] 
whereas using Killing's equation, the symmetric part is 

The terms of ()56p have the following translation into our tensorial language: 
j d'^xTtd^.Li'duL'' = -\ j d'^xU^'^U^'^Kp 

j d^xTTL.U^L^L^ = j d\d,^^d^^^d,^^d^^' {J^T^^^^, + \t'^ 
j d'xTrL.L^L^L'^ = j d'xd.^^d^^^d^^^d^^' (i^^i'^ " ^^£7^ + 

J d^xTV(L^L^)Tr(L^L'^) = d^x d^ip'^d''^p^d^ip^d''ip^T 

J d^xTV(L^L,)Tr(L^L'^) = d^xd^^'^d^'^^d.^^d'^^T^^^ 



(4) 



1) 

1/375 



24 



where Dc/j" = X/^d^ip"^. One can further manipulate the third term integrating by parts and 
commuting covariant derivatives. One finds 



and using (j43p one can further substitute the Riemann tensor by T^'^^ . In the fourth term one can 
ehminate T^^\ 

One has to note that Hasenfratz's action has to be compared to minus our action. This is because 
it appears with the positive sign in the exponent of the functional integral (this is consistent with 
the fact that the (O^p)'^ term has a negative coefficient in (1561) ) . It is then straightforward to 



calculate the following relations between the couplings used in 



12[ | and our couplings: 



5--; A-/; + /2--; h- ; /4 - j ; 

With these relations, one can translate his beta functions and one finds that they agree with 
those given in section HID, with a single exception: the term proportional to a nd containing 
no fi in f3g2 . We observe that the two polynomials in the c's in equation (39) in [12t | are the same, 
up to an overall factor 2. As a consequence, when one extracts the beta function of cq//^ = ^/g^ 
and rewrites it in terms of the /j's, the coefficients of and cp'X are exactly the same. This differs 
from the beta function given in (j53p . where the two coefficients differ in the last term. We believe 
that this difference can be attributed to the different cutoff scheme. 
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